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ABSTRACT 

We investigated the relationship between supermassive black hole (SMBH) mass and 
host starburst luminosity in Seyfert galaxies and Palomar-Green QSOs, focusing on 
the host galaxy morphology. Host starburst luminosity was derived from the 11.3 
[xm polycyclic aromatic hydrocarbon luminosity. We found that the SMBH masses 
of elliptical-dominated host galaxies are more massive than those of disk-dominated 
host galaxies statistically. We also found that the SMBH masses of disk- dominated 
host galaxies seem to be suppressed even under increasing starburst luminosity. These 
findings imply that final SMBH mass is strongly regulated by host galaxy morphology. 
This can be understood by considering the radiation drag model as the SMBH growth 
mechanism, taking into account the radiation efhciency of the host galaxy. 
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1 INTRODUCTION 

Recent observations have found that the mass of super- 
massive black holes (SMBH) in galactic centers correlates 
stron gly w ith bulge mass in both active and inactive galax- 
ies (IKprm cndv fc Gcbhardt 2001; Morritt & Fcrrarese 
200ll: iTremaine el^ al.1 l2002i : iMarconi fc HuntI l2003l : 
Barth. Greene, fc Hoi |2005| ). Also, it has been revealed 
that SMBH mass does not correlate with galaxy disk 
mass; in fact, the SMBH-disk mass r atio is much smaller 
than the SMBH-bulge ma ss ratio (|Salucci et all l2000l : 
iKormendv fc GebhardtlboOll) . These observational findings 
suggest that SMBH formation is strongly connected to the 
bulge component, not the disk component. However, the 
physical mechanism of SMBH formation that can lead to 
these observational results has remained unclear. 

Investigations of starbursts around active galac- 
tic nuclei (AGNs) may prove pertinent to this ques- 
tion. Observations of polycyclic aromatic hydrocarbon 
(PAH) emission have gradually revealed evidence of star- 
burst phenomena around AGNs a nd a connection be- 
tween AGN activity and s t arbur s ts lUmanishil [20021. | 2003| : 
Rodn'Euez-Ardila fc Viesasl '2003"; 'imanishi fc Wada' '2OOJ: 
Schweitzer ct al. 2006; Maiolino ct al. 2007; Nctzcr ct al] 
20071: iShi et al.ll2007l : IWatabe. Kawakatu. fc Imanishill200i : 
Lutz et alT 20081 '). This means that the mass accretion pro- 
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cess onto the SMBH, i.e., the SMBH growth mechanisms, 
could be closely connected with starburst phenomena. 

lUmemural (|200ll ') and iKawakatu fc Umemural (|2002l 'l 
suggested that the physical mechanism of the link between 
the SMBH and bulge formation may be the radiation drag 
effect (the Poynting- Rob ertson effect ) from bulge starbursts. 
Especially, IKawakatu fc Umcmura| (|2004l ) showed that the 
radiation efficiency differs for starbursts in the bulge com- 
pared to the disk. The bulge is round, providing high radi- 
ation efficiency. Disk starbursts are less efficient than bulge 
starbursts in the same starburst luminosity range due to 
photon escape from the disk surface and edge-on opacity 
consideration, so the final SMBH mass of a host galaxy with 
a disk starburst cannot be large. Thus, to understand SMBH 
growth mechanisms, it may be necessary to consider where 
the starbursts occur in their host galaxies. In order to un- 
derstand the radiation effects from the host starburst and 
confirm the radiation drag model, we must investigate the 
relationship between SMBH mass and host starburst activ- 
ity, focusing on host galaxy morphology. 

To date, the morphology of Seyfert galaxies (low 
luminosity AGNs) have been well studied and is almost a ll 
spiral (e.g., iMcLeod fc Riekd [l995l : iHunt fc Malkanlll999l ). 
Moreover, the starbursts in the galaxy disk have been inves- 
tigated (Watabe, Kawakatu, fc Imanishi 2008). For QSOs 
(high luminosity AGNs), although it has been diflicult to 
determine their morphology due to their luminous nuclei, 
recent high-resolution imaging observations in the optical 
and near infrared, and those exploiting adaptive optics 
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Table 1. QSOs properties. Column 1, 5: Object name. Col. 2, 6: 11.3 yim PAH luminosity. Col. 3, 7: SMBH m ass estimated by Kawakatu et 
al. (20 07). t: the data are f rom V ester gaard et a l. (2006). Col. 4, 8: References of host galaxy m orphology; [1] Guvon, Sanders, & Stockton 
1I2OO6I): [2] Iveilleux et al.l 1I2OO6I): [31 iDunlop et al. (20 03); [4] iHamilton. Casertano. fc Turnshek (20 02,); [5] .Marguez et al.. (.2001.) ; [6] 
iMcLeod fc McLeodI 1I2OO1I) ; f7l |Percival et al.l bOOll) ; [8T lBahcall et al.l lll997l) ; [91 iTavlor et al.l ^19m . 
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20031: Veilleux et al.l 20061: ICuvon. Sanders. & StocktonI 
20061). Thev found not only the elliotical (which is eauiva- 



lent to the bulge, classified as the spheroid) component but 
also the prominent disk component in the host of lower lu - 
minosity QSOs and radio-quiet QSOs (|Dunlop et al.ll2003h . 
Also, we now know that QSO host galaxies are gas-rich 
(which is not normal for ellipti cal galaxies) (| Evans et al.l 
I2OOII . I2OO6I : IScoville et al] l2003h and that starbursts oc- 
cur in the host galax i es (iHaas et al.l 20031: iBarthell I2OO6I: 



Schweitzer et al 



20061: iMaiohno et a.1. 



20071 : IShi et al.ll2007l : lLutz et al.ll200^ 



20071 : iNetzer et all 



In this letter, to clarify the SMBH growth mechanism 
that satisfies the observational results, we investigated the 
SMBH mass-host starburst connection, taking into account 
the host galaxy morphology for Seyfert galaxies and QSOs. 
Throughout the paper, we adopted Ho = 80 kms~^Mpc~^, 
Qm = 0.3, and SIa = 0.7. 



2 DATA AND ANALYSIS 
2.1 Sample data 

To investigate the effect of starbursts and host galaxy mor- 
phology on the final SMBH mass, we selected AGNs based 
on estimated SMBH mass, PAH emission, and morphology. 
The samples in this paper are Palomar-Green (PG) QSOs, 



selected at B band to have blue U — B color (Schmidt & 
Green 1983). Since the B- and U-band emission of PG-QSOs 
are dominated by the AGN, there are no biases for the pres- 
ence of starbursts. 

Se yfert galaxies in the CfA llHuchra fc Bur3ll992l ) and 
12 /^m l|Rush. Malkan. fc Spinogliolll993l ) samples, selected 
on the basis of their host galaxy magnitudes and IRAS 12 
Hm fluxes, respectively. These samples do not also include 
some biases for the presence of starbursts. 



2.2 Black hole mass 

To estimate SMBH mass, Mbh, of PG-QSOs and Seyfert 
1 galaxies, we used a method based on the reasonable as- 
sumption that the motion of ionized gas clouds around the 
SMBH is dominated by the gravitational force, and that 
the cl ouds within the broad-line region (BLR) are virialized 
(e.g., iPeterson fc Wandel|[l999l ). The velocity dispersion v 
can be estimated from the FWHM of H/3 broad line emis- 
sion V — /iipwHM by assuming the isotropic spherical virial 
coefficient, / = \/3/2 (Netzer 1990). We selected the H/3 
because H/3 lines radiate by simple photoionization mecha- 
nisms and that their line profiles refiect gravitational poten- 
tial of th e H/3 emiss ion region. Adopting an empirical rela- 
tionship iKaspi et al. 2000) between the size of the BLR and 
the rest-frame optical continuum luminosity, XLx{5100A), 
and using reverberation mapping, we obtain the following 
formula: 



Mb 



: 4.9 X 10 



ALa(5100A) 



1044 ergss-i 



IIPWHM 

,103 kms- 



M0.(1) 
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For the error of SMBH mass, McGill et al. 2008 have com- 
pared 12 formulae taken from the literature, showing that 
SMBH mass estimates can differ on average 0.13 ± 0.05 or 
0.38 ± 0.05 dex in the case of the same or different virial co- 
efficient, re spectively. For Sey fert 2 galaxies, we used SMBH 
mass data (jBian fc Gull2007l ) whose mean error was within 
a factor of 1.6. These were estimate d from the SMBH ma ss- 
stellar velocity dispersion relation jXremaine et al.ll2002h . 



2.3 Host galaxy starburst luminosity 

Since the PAH molecules are excited by far-UV photons in 
the photo-dissociation region around the HIT region, and 
strong PAH emission is often observed from even a weak 
starburst (jlmani shi 2002') , we can use PAH emission as an in- 
dicator of starburst activity. We used the PAH em ission esti- 
mated bv lWatabe. Kawakatu. fc Imanishi^ (l200^ for Sey fert 
galaxies (6.2, 7.7, and 11.3 fj.m) and Shi ct aU (|2007f l for 
PG-QSOs (7.7, and 11.3 fitn). We select the 11.3 PAH 
emission. The 7.7 PAH emission is sometimes affected by 
the broad and strong 9.7 fim silicate absorption, especially 
in Seyfert 2 galaxies. Thus, it could be difficult to distin- 
guish between 7.7 fim PAH emission and 9.7 fim silicate 
absorption. PAH emission was obtained wit h the Spitzer 
Spac e Telescope Infrared Spectrograph (IRS; iHouck et al.l 
120041 ') (| Werner et all 120041 '). For Seyfert galaxies, since the 
PAH emission was observed with the slit-scan mode (PID 
3269, PI: J. Gallimore), the entire host galaxy regions of the 
Seyfert galaxies were covered. For PG-QSOs, PAH emission 
was obtained by the slit width of Short- Low (SL) (SLl: 3". 7, 
SL2: 3".6) and Long-Low fLL) (LLl: 10" 7, LL2: 10". 5) mod- 
ules (for PID and PI, see ishi et alll2007l l. The SL slit width 
is roughly co mparable to the effective radius of PG-QSO 
host galaxies (|Guvon. Sanders, fc Stocktonll200(j ') and sev- 
eral hundred-parsec to kiloparsec-scale starbursts have been 
conside red as the origin of the far-infrared radiation in PG - 
QSOs iHaas et all [20031 : lBarthell2006l : lNet"zer et al.|[2007l ). 
Thus, we consider that these slit observations cover almost 
all starburst activity in PG-QSOs host galaxies. 



2.4 Host galcixy morphology 

We used host galaxy morphology classifications for PG- 
QSOs based on the literature, evaluated by 2D fitting of 
the obtained image (references listed in Table 1). We care- 
fully checked the morphology classifications and defined a 
host galaxy as elliptical- or disk-dominated only in cases 
for which the literature listed said galaxy as only elliptical 
or disk, respectively. We also defined a bulge -I- disk host 
galaxy in cases where a 2D fit favored a two-component 
(bulge + disk) model, or in cases where the morphology 
decision given in the literature varied (see Table 1). We 
also used a homogeneous morphology classification crite- 
rion to check ou r results, by separately examining o bjects 
classified only bv lGuvon. Sanders, fc StocktonI l|2006l ). This 
sample consisted of a number of PG-QSOs (20 objects) in- 
vestigated by near-infrared AO imaging with the Gemini-N 
and Subaru Telescope. For the Seyf ert galaxies, since their 
hosts are almost all s piral galaxies (jMcLeod fc Riek3ll995l : 
iHunt fc Malkanll 19991 ). we assumed that their morphologies 
were disk-dominated. 



3 RESULTS: SMBH MASS AND HOST 

STARBURST GEOMETRY RELATIONSHIP 

We plotted the 11.3 fj.m PAH luminosity and SMBH 
mass specifying the host galaxy morphology in Figure [T] 
The left and the right panels of this figure show that 
the morphology classification is u sed Table 1 and only 
iGuvon. Sanders, fc StocktonI l|2006l ). respectively. 

We applied detailed statistical tests about the difference 
of these distributions (Takeuchi et al. 2009 in prep.) . To 
summarize, the SMBH mass distributions of the elliptical- 
and disk-dominated host galaxies are significantly differ- 
ent. However, the difference of the PAH luminosity is not 
very clear. (We showed this statistical analysis in Ap- 
pendix.) Averaged SMBH mass and its dispersions are 

< logMBH/M© >= 8.42 (±0.44) and 7.48 (±0.36) for the 
elliptical- and disk-dominated host galaxies, respectively. 
Also, bulge ± disk host galaxies are distributed in both the 
elliptical- and the disk-dominated host galaxy regions. In 
particular, for the disk-dominated host galaxy, although the 
PAH luminosity increases by about three orders of mag- 
nitude, SMBH mass increases by only about one order. 
These results indicate that the final SMBH mass is strongly 
connected with host galaxy morphology; the SMBH of a 
disk-dominated host galaxy is suppressed, while a more 
massive SMBH can form in an elliptical-dominated host 
galaxy. In order to remove the differences in observations 
and morphology classification methods, we checked our find- 
ings using data from Guyon, Sanders, fc Stockton (2006|) 
only in the right panel of Figure 1. Although the sample 
size decreases, the tendency of our results does not change; 

< logMBH/Mo 8.40 (±0.41) and 7.42 (±0.34) for the 
elliptical- and disk-dominated host galaxies, respectively. 



4 DISCUSSION: AGN FUELING MECHANISM 

To interpret our findings, we must consider the SMBH 
growth mechanism including both the host starburst effects 
and the host galaxy mor phology. Although gala xy merg- 
ers ( e.g., iHernguist Il989l) and stellar bars (e.g., iNoguchil 
Il988f ) have also been considered as SMBH growth mech- 
anism candidates, the relationship between final SMBH 
mass and these mechanisms is still unknown. Therefore, 
it is difficult to explain the difference in SMBH mass for 
the same starburst luminosity range using these mecha- 
nisms. Thus, in order to relate the host starburst and 
host galaxy morphology with SMBH formation, we focused 
on the radiation-hydrodynamic effect from the host star- 
burst. The radiation drag is a relativistic effect known as 
the Poynting-Robertson effect. It is a possible mechanism 
for extracting angular mo mentum from the gas and driv- 
ings ]\ffiHjnass accretion^ fc Mincshig3 
1 19971: lu memur al l200ll : iKawakatu fc Umemurali2002i ). Final 
SMBH mass is connected with the absorption efficiency of 
the amount of radiation energy from the starburst. This ra- 
diati on drag efficiency is strongly affected by host ge ometry 
(U memura. Fukue. fc Mineshig c 1997 , 1998; Ohsu ga et al.l 



: J I ' ' ' II ' -- i'-' ^^^^^^^ 

[l999|). IKawakatu fc Umemural (|200J) explored the possi- 



bility that SMBH mass of a disk-dominated host galaxy 
could be one to two orders of magnitude smaller than 
that of an elliptical-dominated host galaxy due to ef- 
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Figure 1. The 11.3 /im PAH luminosity vs. S MBH mass for a gi ven the host m orphology. The morphology classification of PG-QSOs 
host galaxy is used Table 1 (left) and only Guyon, Sanders, & Stockto^ ll2006h (right). Open circle, open triangle, and filled square: 
elliptical-dominated, bulge+disk, and disk-dominated PG-QSOs, respectively. Filled diamond: Seyfert galaxies. 



fects of geometrical dilution and opacity. If an elliptical- 
dominated host galaxy begins with starbursts in a highly in- 
homogeneous and optically thick interstellar medium (e.g., 
ISanders et all 1 19881 : iGordon. Calzetti. fc Wittlll997l ). radi- 
ation drag c ould effectively work to extr act the angular 
momentum (|Kawakatu fc Umemura|[2003 ). In contrast, a 
large number of photons escape from the disk surface of 
a disk-dominated host galaxy. Also, radiation from a disk- 
starburst is shielded due to edge-on optical depth. Therefore, 
from a radiation-hydrodynamic point of view, final SMBH 
mass could be strongly connected to starburst location. Also, 
there seems to be no possibility that the radiation efficiency 
is equivalent for disk- and elliptical-dominated host galaxies 
and the SMBH mass of a disk-dominated host galaxy is cur- 
rently small, but will grow up in the future due to SMBH 
growth delay. It's because SMBH masses of the inactive spi- 
ral galaxies are comparable to those of the disk-dominated 
host galaxies with AGNs and follow the same SMBH-bulge 
mass relation (jPastorini et al.ll2007h . Therefore, in the case 
that all inactive spira l galaxies have passed the AGN phase 
l|Marconi et al.l 120041 ') . the final SMBH masses of the disk- 
dominated host galaxies with AGNs could not reach those 
of the elliptical-dominated host galaxies. Thus, our result 
flows naturally from a radiation drag model that includes 
the radiation efficiency due to host galaxy morphology. 



5 SUMMARY 

We investigated the SMBH mass-host starburst connec- 
tion, taking into account the host galaxy morphology for 
Seyfert galaxies and PG-QSOs. We checked the statistical 
test about the difference of these distributions. As a result, 
we found that host galaxy morphology may strongly regu- 
late final SMBH mass, as the SMBH masses of elliptical- 
dominated host galaxies were more massive than those of 
disk-dominated host galaxies. Also, the SMBHs of disk- 
dominated host galaxies showed suppressed mass even in 
the case of increased starburst luminosity. These findings 
indicate that the SMBH growth mechanism is strongly con- 
nected to radiation efficiency dependent on the geometry of 
the host starbursts. 
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APPENDIX 

Here we explained the statistical tests. (Further details will 
be found in Takeuchi et al. (2009 in prep.).) First, we tested 
if the vertical marginal distributions of Figure 1 are the 
same. This is a well-known statistical problem which is 
properly addr essed by the Kolmogorov-Smirnov two sam- 
ple test (e.g., ISiegej1l956l l. The test statistic is the dis- 
tance between the two distribution functions (DFs): D = 
max|_Fa(MBH) — Ff,{MBii)\ , where the subscripts a and b 
represent that they are the DFs of the sample in the left 
and right panel of Figure 1, respectively. Then, it is known 
that a quantity 

^2 _ 4I?^nciiWdisk 
ndi + ndisk 

obeys the x'^-distribution with a degree of freedom (dof) 
of 2, under the null hypothesis that these are drawn from 
the same DF. The distances D from the sample of a and b 
are 0.818 and 0.882. Then, the probabilities having the 
values [eq. corresponding to these values are 5x 10~® and 
0.0001, respectively. Thus, in both cases, the null hypothesis 
is very clearly rejected. 

Next, we examine the horizontal marginal distribu- 
tion of the samples. Along with the abscissa in Fig- 
ure 1, we have a significant number of upper limits both 
in elliptical and disk galaxy samples. In this case, usual 
two sample tests are no longer valid. We first estimate 
DFs of the PAH luminosity, F(Lpa h) of both samples by 
the Kaplan-Meier (KM) est imator (|Kaplan fc Meied HoHsl : 
iKalbfieisch fc Prenticell2002D . The KM estimator is designed 
to obtain the so-called survival function (or survivor func- 
tion), S{t), of a given sample including right censored data 
(i.e., data with lower limits). A survival function is related 
to the DF by the following relation: S(t) = 1 - F{t). 

To construct this estimator, we first mathematically for- 
mulat e the current problem according to lFeigelson fc NelsonI 
j 19851 ). Since most of the astronomical observational data 
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are selected with a certain detection limit, usually a sample 
includes upper limits, instead of lower limits. In statistical 
terminology, lower-limit data are called "left-censored" . Let 
Ti , T2 , • ■ • , denote measurements, where the superscript 
L means "left". If T/' < Ai, where At is its upper limit, 
then we cannot obtain the exact value of Ti but the up- 
per limit Ai. The information available is a combination of 
{Xf ,5f } = {T„ 1} and {A„0} for T,^ > A, and T,^ < A„ 
respectively. Suppose that T/' is a random sample drawn 
from a distribution F^{t) = P{Ti ^ t) where P{E) denotes 
the probability of an event E. If {Ai}i^i^2,...,n are mutu- 
ally independent, identically distributed (referred to as IID) 
and independent of the true measurements {r/'}i=i,2,...,n, 
this statistical model is called random censorship. We can 
easily transform left-censored data to right-censored ones 
by setting a constant M; Ti = M - Tt , Xi = M - Xt, 
Ci = M - Af, Si = 5f . Then, {Xi,Si}i=i,2,...,n represent 
right-censored data. The KM estimator has the following 
form: 



1 1 i\x' 



t > X 



(3) 



where Xq-^ denote distinct, ordered observed values in 
which ties are identified, and rij = # |A:; ^ a;'y-)}, 

Xk = } . A detailed but corn p rehen sive 
derivation is found in , e.g., If^igelson fc NelsonI l|l985h or 
iKalbfleisch fc Prentice! ||2002D . 

Once we get the KM estimate of the survival function, 
S{t), we can convert it into the estimate of the DF, F{t), as 

= P [M -T ^t) 
= P{T ^ A4 -t) 

= S{M-t). (4) 



Then, 

F^{t) = S{M -t) 



(5) 



In this analysis, we set M = 0, i.e., we made a flip of a sign 
t ~t. 

In order to see if these two DFs are different, we should 
perform a statistical test. Though there are a few methods 
to do so, we adopt the Mantel-Haenszel (MH) two-sample 
test (also often referred to as the logrank test, but the lat- 
ter denotes some variation s ). Th e derivation is shown in, 
e.g.. IKalbfleisch fc Prentice! (|2002! ). The MH test makes use 
of the values as follows: first we sort all the galaxy sample 
along with the (minus logarithmic) luminosity, — log Lpahj- 
Then, we define mean numbers of galaxies at each luminos- 
ity — logI/pAH,j for both samples, rrikj = dj^^^^ where k 
denotes the label of the two samples, "ell" and "disk", and 
rij = ndi.j + ?idiskj-We also define variance 



I'kk,! — 



(fc and 



nkjntjdj {rij — dj) 
(njf {rij - 1) 

£ again denote 



(6) 



"eU" 



and 



"disk" 



and covarianc^; Vke,j = —Vkk,j- Consider 
a deviation vector Aoiij, Adiskj^ = 

doll J — "loll J, ddisk.j — "idiskj,^ • Since the dof is 
2 — 1 = 1 in the current problem, Aoiij = — Adisk.j- 
Then, simply we can use a statistic Ak^/vkk- If the two 
samples are drawn from the same DF (this is the null 
hypothesis to be tested), this statistic should obey the 
X^-distribution with dof = 1. From the sample of a and 6, 
we have ~ 0.694 and 0.951, respectively. Probabilities 
of having these values under the x^-distribution are 0.405 
and 0.329. Then, if we set the confidence limit of 0.05, we 
cannot reject the null hypothesis, i.e., the two DFs are not 
significantly different. 
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